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^ ■ Abstract 

■ The stability of hot expanded nuclear droplets against small bulk and 

| surface oscillations is examined and possible consequences for multifrag- 

mentation are discussed. 

> : 

^ ; 1 Introduction 

O 

Spinodal instabilities were suggested by Bertsch, Siemens and Cugnon ]]] al- 
ready 15 years ago as a possible mechanism leading to multifragmentation of 
Q\ . hot expanding nuclear matter. 

! In heavy-ion collisions we expect that hot nuclear droplets are formed, 

! which subsequently expand leading to low densities in the interior. Below 

certain values of the density, bulk and surface instabilities may occur and lead 
to multifragmentation. 

In this contribution we report on a study of bulk and surface instabilities 
• i-H ', of spherical nuclear droplets as function of temperature and density. 

X ' 

.F*. 2 Bulk instabilities 

We consider the normal modes of a nuclear droplet. Let us define a complete or- 
thonormal set of functions ^\(r, (p) = Afi n ji(k n ir) y™^, <p), where yp(^, <p) 
are modified (real) spherical harmonics, ji(kr) are spherical Bessel functions 
and Min are normalization constants. As we aim to consider the distortions of 
a spherical droplet with radius Ro, we impose the condition ji(k n iR Q ) = 0. We 
consider the distortions defined by the irrotational displacement field 

s(r, t)=Vj2 9a(*) ^A(r) = V w(r, t) (1) 
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Table 1: Parameters of SkM* and SIII Skyrme forces 



to t 3 h t 2 xq x 3 a 

SkM* -2645 15595 410 -135 0.09 1/6 
SIII -1128.75 14000 395 -95 0.45 1 1 



from the general form of the dispacement potential w(r,t) = J2x <l\(t) ^\(r) . 
Note that the surface is not kept fixed with this definition of the displacement 
field. The density varies according to the continuity equation 

^^- + dw[ 6 (f,t)v(r,t)] = 0, (2) 

which guarantees exact conservation of mass, and - together with the condition 
^/\(Ro, Q) = - also of the center of mass. 

In harmonic approximation small oscillations around equilibrium are de- 
termined by the set 

E s av<7a' + £C av <?v = o (3) 

A A 

of coupled equations. The eigenmodes are obtained from diagonalizing 

(C X x> ~ uj 2 B xx ,)q x , = . (4) 

For u 2 > the corresponding mode is stable (q x ~ sin(u;£)). Exponential 
instability (q x ~ exp(yt)) occurs for oj 1 = — -y 2 < 0. 

Analytic expressions for the mass (B xx >) and stiffness (C X y) tensors are 
derived from the velocity field s and the Skyrme energy-density functional , 
respectively. The mass tensor is diagonal. The stiffness tensor C is the sum 
of the contributions C v , C T , C w , C s and C c resulting, respectively, from the 
volume, intrinsic kinetic energy, Weizsacker, surface and Coulomb terms in the 
energy density. The analytic expressions of the tensors are given in JJ. Modes 
belonging to different I, m (multipoles and their components) are decoupled. 
The only couplings left are those corresponding to different numbers n of nodes 
in the displacement field for the same multipolarity. The eigenvalues are ob- 
tained by numerical diagonalization of the matrix B^ 1 C. The calculations 
have been performed for two Skyrme forces, i.e. SkM* and SIII implying, re- 
spectively, a soft and a stiff equation of state (EOS). The parameters of both 
forces are listed in Table |]. The contribution from the intrinsic kinetic energy 
has been calculated in the adiabatic limit (constant entropy with isotropic 
momentum distribution). 
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Fig. 1. Borders of spinodal instabilities for infinite symmetric nuclear matter 
(heavy solid line), for a gold nucleus (A=197, dashed line) and for a smaller 
nucleus (N=Z=50, dotted line) as calculated for the soft (left) and the stiff 
(right) EOS, respectively. Temperatures T are given in MeV. 

From now on, for all figures, p/po denotes the ratio of the actual density 
to normal nuclear density (p = 0.16 fm~ 3 ). Fig. 1 displays the areas of 
spinodal instability for the density modes with / = 2 for three cases, i.e. infinite 
symmetric nuclear matter, a gold nucleus (Z=79, A=197), and a symmetric, 
roughly two times smaller fragment (Z=N=50) calculated with SkM* (left) 
and SIII (right) forces, respectively. For infinite nuclear matter the area of 
spinodal instability is substantially larger than those for finite systems. The 
difference between the spinodal lines of the two finite systems is small. The 
areas of spinodal instability for the higher multipoles 1 = 3,4,... are further 
reduced. The results for infinite symmetric nuclear matter have been obtained 
by neglecting surface, Coulomb and Weizsacker terms (Cj^) = C^J = C A ^P = 
0) and taking the limit (R — > oo) in the volume and kinetic-energy terms 

l°AA' > °AA'J- 

Fig. 2 shows quantitatively the importance of different contributions to the 
lowest eigenfrequency as function of the density at a typical temperature of 
4 MeV. Due to their small values, the surface and Coulomb terms have prac- 
tically no influence on the stability of the nuclear droplet. There is a delicate 
balance between the contribution from the kinetic energy term C^y and the 
volume term CyQ, such that the role of the Weizsacker term Cj^JP becomes 
crucial. This term grows substantially with density and thus is important for 



the increasing stability with increasing density. 




Fig. 2. Contributions of different terms to the lowest diagonal element of the 
B _1 C matrix at T = 4 MeV for a gold nucleus and soft (left) and stiff (right) 
EOS. Additionally the lowest eigenvalue to 2 is displayed by the solid line. The 
difference of this eigenvalue from C tot /B is due to couplings to higher-n modes. 
Note that contributions from Coulomb and surface terms are very small. 

The crucial quantity in the multifragmentation process is the instability 
growth rate 7 (for u 2 < 0, q\ ~ exp(7t)) or the characteristic growth time r = 
I/7 for a particular mode. Multifragmentation, initiated by such instabilities, 
can occur only if these characteristic times are short enough compared to the 
characteristic evolution time of the system. In Fig. 3 (left) we present the 
shortest characteristic times as function of density and temperature for the 
bulk instabilities of the gold system calculated with the soft EOS. 

3 Surface instabilities 

For early stages of the expansion, where densities and temperatures are still 
high, no bulk instabilities exist. There instead, e.g. for initial temperatures 
higher than 8 MeV, the surface vibrations are found unstable. However, the 
characteristic times for these instabilities are about an order of magnitude 
larger than those for bulk oscillations in accordance with @]. We follow here 



the standard Bohr and Mottelson theory ||. Fig. 3 (right) presents the small- 
est characteristic times for these surface instabilities. Although our basis is 
complete and - in principle - a linear combination of the collective displace- 
ment fields ty\ can describe an arbitrary collective motion, including surface 
vibrations of incompressible matter, such surface motion requires a very large 
number of terms due to slow convergence. Therefore, we calculated the surface 
instabilities separately. As we see from Fig. 3 bulk instabilities are dominant 
at small densities and temperatures, while for large enough temperatures and 
densities only surface modes are unstable. 




Fig. 3. Minimal characteristic times for bulk instabilities of a gold nucleus 
calculated with a soft EOS (left 3d-plot). All modes with I = 2,3,4,5 and 
n = 1, . . . 6 are taken into account. The contour lines correspond to the values 
20, 25, 30, 40, 60 and 100 fm/c. The temperature T is given in MeV. The 
right 3d-plot shows the same (rotated by n around the vertical axis) for surface 
instabilities. Here, the contour lines correspond to the values 150, 200, 250, 
300, 400 and 600 fm/c. Note that in both diagrams the system is stable in 
the regions outside the holes. 

4 Consequences for multifragmentation 

The expansion of hot nuclei has been studied in refs. || for soft and stiff equa- 
tions of state and compared with experimental data J7|. For the soft EOS the 
expansion trajectories in the (T, p) plane reach turning points around T 4 
MeV and p/po ~ (0.25 — 0.45). Around these turning points the collective 
motion is very slow, and hence fast enough unstable modes can develop and 
initialize multifragmentation. Although finite size effects substantially reduce 
the area of spinodal instabilities with respect to that of infinite nuclear matter, 



it is clear from Fig. 1 that the turning points essentially remain in the insta- 
bility regime. From Fig. 3 (left) we see that characteristic times for density 
instabilities around this turning points are on the order of 25-40 fm/c. Such 
times seem to be short enough to cause an irreversible decay of the system 
towards multifragmentation. Usually several modes become unstable in this 
region, with close characteristic times, such that the production of many dif- 
ferent fragments is possible. Our study of surface modes show that initially 
deformed nuclear droplets can hardly become spherical, because the charac- 
teristic times for shape restoration is large compared to the expansion time. 
For initial temperatures larger than 8 MeV such initial deformation will tend 
even to grow. 
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